Numerical solutions of incompressible Navier–Stokes equations using modi ed Bernoulli’s law

Simulations of incompressible ows are important for many practical applications in aeronautics and beyond, particularly in the high Reynolds number regime. The present formulation is based on Helmholtz velocity decomposition where the velocity is presented as the gradient of a potential plus a rotational component. Substituting in the continuity equation yields a Poisson equation for the potential which is solved with a zero normal derivative at solid surfaces. The momentum equation is used to update the rotational component with no slip=no penetration surface boundary conditions. The pressure is related to the potential function through a special relation which is a generalization of Bernoulli’s law, with a viscous term included. Results of calculations for twoand three-dimensional problems prove that the present formulation is a valid approach, with some possible bene ts compared to existing methods. Copyright ? 2003 John Wiley & Sons, Ltd.